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The intersection between dislocations and a Ag(lll) surface has been studied using an interplay 
of scanning tunneling microscopy (STM) and molecular dynamics (MD). Whereas the STM provides 
atomically resolved information about the surface structure and Burgers vectors of the dislocations, 
the simulations can be used to determine dislocation structure and orientation in the near-surface 
region. In a similar way, the sub-surface structure of other extended defects can be studied. The 
simulations show dislocations to reorient the partials in the surface region leading to an increased 
splitting width at the surface, in agreement with the STM observations. Implications for surface- 
induced cross slip are discussed. 
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Dislocations at surfaces play a major role in many ar- 
eas of materials science. Dislocations intersecting sur- 
faces are for example very important for the control of 
growth and solidification processes. Steps terminating at 
screw dislocations are ideal for continued crystal growth 
and dislocations can strongly modify the surface stress 
important for atomic mobility and island nucleation in 
growth from the gas phase ||^ . Dislocation-surface inter- 
actions also play a crucial role in determining the fracture 
toughness of a material since the emission and absorption 
of dislocations at a crack tip controls the possible blunt- 
ing of the tip ||]. 

However, the atomistic understanding of the behavior 
of dislocations at surfaces is still very scarce ^, || |^, 

|9|. Bulk dislocations have in the past been studied 
extensively with electron microscopy (EM) and in some 
cases it has been possible to obtain atomic-scale informa- 
tion by imaging columns of atoms along a straight bulk 
dislocation It is not possible to obtain similar reso- 
lution with EM for dislocations at surfaces. STM on the 
other hand is the technique of choice to reveal atomic- 
scale surface structures such as surface dislocations, but 
STM is blind with respect to the region beneath the top- 
most layer. 

In this Letter we show how STM observations of a 
dislocation intersecting a surface can be combined with 
atomistic simulations to provide detailed atomistic in- 
formation about the structure of the dislocation in the 
near-surface region. The good agreement between ex- 
perimentally observed and simulated surface structures 
gives credibility both to the interpretation of the surface 
structures as signatures of dislocations and to the calcu- 
lated structures below the surface. The same principles 
can be used to study the interactions between surfaces 
and other extended defects such as grain boundaries. 



Consider a bulk dislocation that intersects a surface. If 
the total Burgers vector has a nonzero component along 
the surface normal, a surface step will end at the disloca- 
tion. If the dislocation is split into partial dislocations, a 
step may be seen between the partial dislocations, even if 
the total Burgers vector is parallel to the surface. From 
an atomic resolved STM image, the in-plane component 
of the Burgers vector can also be determined, allowing a 
direct determination of the Burgers vectors of the indi- 
vidual partial dislocations. 

The STM cannot, however, determine the edge or 
screw character of the dislocations as it gives no informa- 
tion about the line vector of the dislocation. Information 
about the line vector and the sub-surface structure can 
be obtained from an interplay with MD simulations. 

The STM measurements were performed on a single 
crystal Ag(lll) surface in two ultra-high vacuum sys- 
tems equipped with standard facilities for sample prepa- 
ration and characterization. One chamber houses a fast- 
scanning, variable-temperature STM the other one 
a low-temperature STM ||l2|. The clean Ag(lll) surface 
was prepared by several sputtering-annealing cycles. 

Several kinds of dislocations were found in the studied 
samples. In a number of surface step ends at 

a dislocation, see Fig. |l|. Despite a significant amount 
of surface diffusion, the structure of the dislocation itself 
is very stable. Up to a step height of about 0.13 nm, 
i.e., around 2/3 of the full step height, the positions of 
the atoms do not change within half an hour. Above 
this step height the step shows the frizziness expected for 
Ag(lll) at this temperature, as the step is moving while 
the STM data are collected |||. This inhibits the 
determination of the step profile with the same precision 
as on the lower parts of the step. 

In rare cases dislocations with a Burgers vector in the 
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FIG. 1: Atomic resolution image of the intersection of a bulk 
dislocation with the Ag(lll) surface and its surroundings (in- 
set). The Burgers vector is ^(110). For better visibility, the 
contrast is enhanced on a broad stripe in the middle of the 
image to both the left side and the right side of the step, thus 
no tip change occurs. The horizontal shift in the atomic po- 
sitions on the two sides of the extended dislocation is due to 
the in-surface component of the Burgers vectors of the par- 
tial dislocations. The STM parameters were U = —0.22 V, 
/ = 1.4 uA, and T = 330 K. 



surface plane were observed, and in these cases only a 
minor perturbation of the surface is seen. In a few cases 
some features which we interpret as Lomer-Cottrell locks 
(sessile dislocations with edge character were seen. 

The simulations were performed using a previously 
determined Effective Medium Theory (EMT) potential 
ps] , ^ for Ag. The potential has been fitted to the 
elastic constants, the vacancy formation energy, and the 
intrinsic stacking fault energy in accordance with exper- 
imental and ab initio data JTtI . The simulations have 
been done at zero Kelvin by minimizing the energy of 
the system, and at 300 K using Langevin dynamics with 
a time step of 5 fs . All simulation cells are rectan- 
gular with (111) surfaces at the top and bottom of the 
cells, with free boundary conditions on all surfaces. It 
is not possible to have periodic boundary conditions in 
the directions perpendicular to the surface due to the net 
Burgers vector in the simulation cell. Care was taken to 
insure that the systems were large enough to prevent un- 
wanted interactions with the other surfaces of the system. 

Fig. ^ shows two calculated equilibrium configurations 
of a dislocation with a Burgers vector of 6 = ^[110]. 
Fig. ^(a) corresponds to the lowest energy. The disloca- 
tion can dissociate according to one of the reactions 




FIG. 2: Two configurations of a dislocation with b = 
i[110], corresponding to two different local minima of the 
energy. The simulated system contains approximately 2 mil- 
lion atoms. The atoms in local fee order and the atoms on 
the surfaces perpendicular to the (111) surfaces of the simu- 
lation cell have been removed, corresponding to showing 1 % 
of the atoms. The atoms in local hep order are dark while the 
atoms in the dislocation cores are grey. Configuration (a) has 
the lowest energy. Configuration (b) has a constriction where 
the slip plane of the dislocation changes, and this is shown in 
the blow up. Notice how the separation between the partials 
changes near the surfaces. 



110] ^ i[121] + i[211] dissociating on (111). (2) 



110] ^ i[211] + i[121] dissociating on (111) (1) 



The simulation in Fig. |^(a) was set up as two partials 
according to reaction ^ with the partials initially sepa- 
rated by 2 nm at the surfaces, and the configuration was 
then allowed to relax. The equilibrium configuration has 
a line vector of [110] and is thus a screw dislocation, al- 
though this does not give the shortest dislocation length. 
The dislocation of Fig. ||(a) was initially set up with this 
line vector, but dislocations with other line vectors were 
seen to rotate to this orientation, starting at the surfaces. 

Linear elasticity theory predicts that a screw disloca- 
tion has lower energy than an edge dislocation and that 
the elastic repulsion between two screw dislocations at 
right angles vanishes ||l|. The elastic energy is thus sig- 
nificantly lowered if the partials rotate away from the 
bulk orientation to obtain a more screw-like character. 
This process is only possible because the surfaces break 
translation symmetry along the dislocation line, result- 
ing in the changes seen near the surfaces in Fig. ||(a). 
Similar effects have been seen in simulations of a screw 
dislocation intersecting a Cu(llO) surface [ p^ . 

The surface imprint of the dislocation at the upper sur- 
face in Fig. ||(a), where the distance between the partials 
is increased, corresponds well with the observed struc- 
ture in Fig. |l|. The profile of the surface step ending at 
the dislocation is shown in Fig. |^ for both the experiment 
and the simulation, and the splitting width as well as the 
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FIG. 3: Step profiles extracted from the STM image in Fig. ^ 
(soUd line), from the simulation in Fig. ^(a) {T — OK, dotted 
line), and from a similar simulation at T = 300 K (dashed 
line). 



widths of the individual partials agree well. 

The direction of the surface step is not determined by 
the dislocation, and it may change if surface diffusion 
adds or removes atoms at the step [Q. In the simula- 
tions, the surface step was placed in the same way as in 
Fig. |l|, i.e. adjacent to the partial with the largest out-of- 
surface component of the Burgers vector. If an extra half 
layer of atoms was added to the surface, the step would 
be in the opposite direction and we would see a high step 
followed by a low step in Fig. |[ 

Another sub-surface structure is possible for the dis- 
location, that of Fig. ^(b), obtained when the energy 
of the perfect (undissociated) screw dislocation is mini- 
mized. In this case, dissociation begins at the surfaces 
according to reaction at the top surface and reaction 
H at the bottom surface, in both cases creating partial 
dislocations rotated towards screw character. The screw- 
like constriction formed in Fig. ^b) has a low energy. In 
Cu such a constriction has been shown to have negative 
energy compared to the straight dislocation p9| . 

The energy difference between the two configurations 
in Fig. H is 4.5 eV, favoring the unconstricted configura- 
tion in Fig. H(a). The surface imprint of the constricted 
dislocation in Fig. ||(b) is not seen in the STM. Notice, 
however, that TEM studies of dislocations in thin foils 
of Cu 10 at.% Al ^ revealed both structures in 
Fig. ^ Also, Fig. ||(a) indicates that a second much lower 
splitting width of about 1 nm with no discernible plateau 
between the partials should be observed experimentally, 
corresponding to the configuration at the lower surface 
of Fig. ^(a). If the dislocation had been dissociated on 
the (111) plane instead of the (111) plane, the narrow 
and wide ends would have been swapped and both con- 
figurations are therefore expected at the same surface. 
We do, however, not observe the narrow configuration, 
except for dislocations with Burgers vector in the surface 



plane (see Fig. where both configurations are seen. 

The altered dislocation structure near the surface may 
influence cross slip rates of screw dislocations. Cross slip 
is the process by which a screw dislocation changes glide 
plane and is an important process in metals in the late 
stages of work hardening . As the cross slip process 
requires two partial dislocations to first recombine into 
a perfect dislocation, the rate depends critically on the 
separation of the partials. The dislocation in figure ||(a) 
will have a much increased probablility of cross slip near 
the bottom surface, but not near the top. The dislocation 
in figure ||(b) has effectively already begun the cross slip 
process, which can proceed by moving the constriction 
down along the dislocation. 

Figs, ^(a) and (b) show dislocations with Burgers vec- 
tors in the surface plane at two different temperatures. 
No surface step ends at the dislocation, but it is never- 
theless visible as it splits into partials that each have a 
Burgers vector with a component perpendicular to the 
surface. The dislocation in Fig. ^(b) crosses a step, but 
this is of minor consequence since the elastic field of a 
step is small. 

In Fig. ^(d) we show step profiles along the disloca- 
tions in Figs. |^(a) and (b) as well as from simulations 
aX T — 300 K and T = K. The dislocation dissociates 
according to the reaction ^[110] i[12l] -h i[2il]. Even 
though the height of the experimental curve could not 
be determined, good agreement is found between simula- 
tions and experiment. The agreement is less good at the 
higher temperatures. Part of the experimental broad- 
ening is due to lower resolution of the image itself and 
especially of the tip. While the step edge in Fig. ^(b) 
has an apparent width of 1 nm (consistent with a perfect 
tip), the step edge in Fig. ||(a) has an apparent width of 
2 nm. 

The splitting width of the dislocations, i.e. the distance 
between the centers of the two partial dislocations, can 
easily be obtained from the step profiles. For the dislo- 
cation in Fig. we obtain a splitting with of 6.4 nm, to 
be compared with 5.5 nm and 7.2 nm in the simulations 
at K and 300 K, respectively. The experimental data is 
insufficient to reliably extract splitting width for the dis- 
locations in Fig. |[ the splitting widths in the simulations 
are 8.1 nm and 11 nm at K and 300 K, respectively. 
All these values are significantly higher than the bulk 
splitting widths, measured to be ds — 2.1 nm for screw 
dislocations in silver, and ~ 8.5 nm for edge disloca- 
tions |2j]. In our simulations, the splitting widths far 
from the surfaces are ds — 1.5 nm and dc — 6.6 nm, in 
reasonable agreement with the experimental values. 

Two mobile dislocations may collide and form a Lomer- 
Cottrell (LC) lock which is a sessile edge dislocation that 
splits on two different planes. This can, e.g., occur ac- 
cording to the reaction 



[101] + i[011 



110]^ 1 [112] -fi [110] + i [112] 
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FIG. 4: Observations and simulations of dislocations and ob- 
servation of a Lomer-Cottrell lock. In (a) and (b) we show 
STM images of dislocations with Burgers vectors in the sur- 
face plane, emerging at the arrows. The temperature was 
T — 318 K and T = 7.5 K, respectively. The semi-circle in 
(a) is a vacancy island. The dislocation in (b) crosses a sur- 
face step. The image constrast has been chosen such that 
greyscales on the two terraces are identical. The simulated 
sub-surface structure of such dislocations at T = K is shown 
in (c). The step profiles of the dislocations (thick lines, (a) at 
the top, (b) at the bottom) and the simulated step profiles at 
T = 300 K and T = K (thin lines, top and bottom, respec- 
tively) are shown in (d). The upper curves are shifted up- 
wards, and the low-temperature experimental curve has been 
shifted to lie on the simulated curve, since the total height 
could not be determined in the experiment. In (e) we show 
a structure consistent with a Lomer-Cottrell lock pinning a 
surface step. 

with the line vector along [110]. A surface structure con- 
sistent with such a dislocation pinning a surface step is 
seen in Fig. ^(e) . The angle of the sharp v-shape is consis- 
tent with the angle between the splitting planes of a LC 
lock. The step height changes by one third at the two 
Shockley partials, since they have a component of the 
Burgers vector perpendicular to the surface. Although 
there is complete symmetry between the two Shockley 
partials, the lowest energy configuration spontaneously 
breaks this symmetry. Linear elasticity theory predicts a 
ratio between the splitting distances of 3.8:1 |^. A lower 
ratio is clearly seen in Fig. ^(e), probably from the in- 
fluence of the surface. Unfortunately, the large splitting 



distances make an atomistic simulation with a realistic 
potential prohibitively expensive. 

In summary, we have recorded STM images of dislo- 
cations intersecting the Ag(lll) surface and measured 
step profiles of the resulting surface steps. These STM 
experiments have been compared to atomic-scale simu- 
lations with good agreement, even on the quantitative 
level. We find that this combination of STM and atomic- 
scale simulations provide a powerful method for studying 
the surface-induced structural changes of crystal defects. 
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